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THE SOLAR WIND 

by 

J. R. Herring and A. L. Licht 
SUMMARY 

Parker f s model of a spherically expanding corona, the "solar wind," is compared 
with D. E. Blackwell’s observations of the 1954 minimum equatorial corona. A significant 
discrepancy is found between the predicted and the observed electron densities at dis- 
tances from the sun greater than 20 solar radii. Blackwell’s data are found to be con- 
sistent with a model in which the corona expands mostly within a disk less than 25 solar 
radii thick, lying within the sun’s equatorial plane. The thickness of the disk as a func- 
tion of distance from the sun is qualitatively explained in terms of magnetic pressure. 

The solar wind is found to have a considerable effect on the lunar atmosphere. First, 
the calculated density of the lunar atmosphere is greatly reduced by collisions with pro- 
tons in the solar wind. If the flux of particles in this wind has the conventional values 
ranging between 10 10 to 10 11 cm' 2 sec' 1 , the calculations yield a lunar pressure of 10" 13 
atmosphere of argon, in agreement with the value predicted by Elsmore and Whitfield on 
the basis of observations on the occultation of radio stars. Second, following a suggestion 
by Gold, it was found that the collisions of solar-wind protons with the lunar surface pro- 
duce an atmosphere of cold neutral hydrogen with a density of 10 5 /cm 3 at the lunar sur- 
face. The density falls off at greater distances in accordance with the inverse- square law 

Estimates indicate that the interaction of solar particles with the neutral hydrogen 
will produce an extended lunar ionosphere with a density of the order of 400 protons/cm 3 
in the vicinity of the moon. 
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THE SOLAR WIND * 


IN 


PROPERTIES OF THE SOLAR WIND DURING 
SUNSPOT MINIMUM 

Recently Parker (Reference 1) suggested that the solar corona may be in a state of 
continual expansion, which he refers to as the fT solar wind.” The Parker model appears 
to effect a reasonable description of the corona’s properties during sunspot maxima. In 
this section the applicability of the Parker model during a sunspot minimum is investi- 
gated. The electron density measurements obtained by Blackwell from the 1954 solar 
eclipse (Reference 2) and from airplane observations of the inner zodiacal light (Refer- 
ence 3) are used. 


Equations of Motion of the Corona 

In the presence of a magnetic field B and a gravitational potential V the steady-state 
equations of motion of an ionized fluid with velocity u, mass density p, and pressure P 
are: 

v • P U = 0 (1) 

expressing the conservation of mass, and 

P U • Vu = - Vp - pVv + ~ V X B X B (2) 

expressing the conservation of momentum. There are also the equations of flux 
preservation 


V X (u X B) = 0 


( 3 ) 


and flux continuity 


V • B = 0 , 


and the equation of state 


P = N k T , (4) 

in which T is the temperature and N is the particle density. The potential v is 

^Originally presented at the First International Space Symposium sponsored by the Com- 
mittee on Space Research (COSPAR), Nice, France, January I960 
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V=-^, (5) 

r 

where G is the universal gravitational constant, m® th e sun's mass, and r is the radial 
distance to the sun. 

The densities p and n can be related to the electron density n if it is supposed that 
the corona consists of one-fifth doubly ionized helium and four-fifths ionized hydrogen. 


p - a M h n 

N = 2jfln, 


(6) 


where m h is the mass of a proton and <*,/3 are the constants 

a = 1.286 , 
f3 = 0.929 . 


The Parker Corona Model 

Parker solved a simplified form of the foregoing equations and derived the velocity 
distributions for a series of assumed temperatures. His calculations have been repeated, 
and the density distribution corresponding to the observed temperature has been derived. 
The temperature also has been calculated independently by fitting his results to the 
observed electron density near the sun's surface. 

In the Parker model a spherical expansion is assumed at a constant temperature, 
and the magnetic field is neglected. With these assumptions, Equation 1 can be integrated 
immediately to give 

n r 2 u = constant ; CO 


and Equation 2 becomes 


aM H Bu 2 
2 " Br 




GM 0 iM H n 


The abbreviations 


m h u 2 


* = T T 


kT ’ 


( 8 ) 

( 9 ) 


X 


r 



where 


n = n o v . 


r Q = the sun’s radius, 


= r 


o » 


n = n at r 
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permit Equation 7 to be written as 

X 2 V \Tq = 

and with this, Equation 8 becomes 

4 -JL 

d£ _ &_ *. (11) 

dx ~ x - 1 

A sketch of the </>- x plane is shown in Figure 1. It can be seen that there is a saddle 
point at <p = 1, x = b/4. Parker picks the solution that goes through the saddle point from 
low to high velocities. This is presumably the only stable solution. The other solutions 
require either too high velocities at r = r 0 , or too high densities at r =<*>. This fixes the 
constant of integration, and from Equation 5, 

<P - In 0 = A - 3 + 4 1n-£ . (12) 

With Equation 10 this is sufficient to determine the velocity and density, once b is known. 

When x is near unity, u is very small. It is then possible to determine b from the 
observed density curve by using the hydrostatic equation valid for zero velocity. 

v = exp [f(i- x )_ ■ (13) 

From Equation 13 and the values for v given by the Jager (Reference 4) for x < 3 , the 
value b = 18.70, which corresponds to a temperature of T = 1.70 * 10 6 °K, is obtained. 

Comparison of the Parker Model with Observation 

The solid curve in Figure 2 shows the velocity calculated on the basis of the Parker 
model. It rises from 1.25 km/sec at the base of the corona to 560 km/sec at the orbit 
of the earth. 


Figure 1 - The 4>-x plane 





4 



Figure 2 - The velocity distributions 


The solid curve in Figure 3 shows the electron density variation predicted by the 
Parker model. The dashed curve shows the density variation as observed by Blackwell. 
The two dotted lines show the margins of error estimated by Blackwell as the result of 
experimental error and the uncertain polarization of the f corona. 

At 150 solar radii the Parker model predicts one-thirtieth of the observed density. 

Nonspherical Model 

It is possible to reconcile Blackwell's results with the concept of an expanding corona 
if the hypothesis of a spherical expansion is dropped. From eclipse observations, it is 
well known that the minimum corona is greatly extendee in the equatorial plane. It will 
therefore be assumed that most of the expanding gas mcves nearly parallel to this plane, 
within a disk of thickness g (g is normalized so that g '1) = l). 

If u is interpreted as the mean radial velocity within this disk, then Equation 8 can 
still be used. The first integral of this equation is the Bernoulli relation 

2 ~ u " 2k/3T (■£ " ln = — v ° ' k/3Tb • (14) 

Equation 7, however, must now be replaced by 


0 * 


X g(x) \1V ~ U { 


(15) 
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The subscript "0" refers to conditions at the base of the corona. 

These equations are sufficient to determine u and g as functions of x if Blackwell’s 
values for v are used and a value for u a is assumed. The choice of u „ is based on the 
condition that the corona expands spherically near the sun. On this basis Uo has the 
value given by the Parker model. 

The velocity u is plotted as the dashed line in Figure 2. It is seen to be nearly con- 
stant at 380 km/sec near the earth’s orbit. The thickness g is plotted as the solid line 

m Figure 4. The dotted lines in Figures 2 and 4 show the margins of error to be expected 
from the error in 

The thickness of the disk rises linearly to about x = 10, reaches a peak at x = 25, and 
falls slowly thereafter. If initially the disk is 1 solar radius in thickness, then near the 
earth it will have a width of 5 radii. 

Magnetic Pressure 

It is possible to obtain a qualitative explanation of the contraction of the disk in terms 
of magnetic pressure. Suppose that the magnetic field B within the disk is due to the 
blowing away of a dipole field. In order to obtain the order of magnitude of B, assume 
that it is nearly parallel to u. By writing 


B = K n u 

it is found that flux continuity requires that K be constant along a streamline, 
using Equation 15, it is possible to write 


(16) 

Thus, by 


xg ‘ (17) 

Here B 0 is a function of the streamline and is determined by the original dipole field at 
the base of the corona. For instance, a field of 10" 5 gauss at x = 150 would require a 
value of B 0 at the sun's equator of only 10 2 gauss. 


These field strengths are sufficiently low to permit neglect of the gradient of the 
magnetic pressure b 2 /8t 7 along the streamlines in comparison with the gradient of kinetic 
pressure. The inclusion of a magnetic field of this magnitude will not seriously affect the 
validity of the Bernoulli Equation (Equation 14). However, the magnetic field can produce 
a force transverse to the streamlines, which is assumed to be balanced by a small trans- 
verse pressure gradient; that is, 


V 


8tt 


V.P. 


(18) 


If Equation 18 is integrated between the center c and the edge e of the disk, holding x 
fixed, then 


1 (B/ - B r 2 ) 


8tt 


p - p 



or 


(19) 
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From Equation 17, 


AB 2 

877 


AB„ 


S77(xg) 2 


( 20 ) 


It can also be assumed that 


Thus 


AP ^ APp 

P - P 0 • 


AB„ 


AP„ 


877(xg)' 


p = AP 0 V , 


AB o' 

877 


= AP a 


and ap q may be eliminated* Solving for g gives 


g - 1 

X w 


( 21 ) 


( 22 ) 


(23) 


This curve is plotted as the dashed line in Figure 4; it has a peak at x = 15 in quali- 
tative agreement with the curve of thickness calculated from the density distribution, but 
it is too high for x < 10. 


Discussion 


Nothing has been said thus far about the polar corona because of the lack of reliable 
density measurements at large polar distances from the sun. It is possible that the polar 
corona is in a state of spherical expansion; in that case its density at large distances will 
be considerably less than in the equatorial corona. 

The sun's equator is known to be inclined at an angle of 7°15 ! to the plane of the 
ecliptic. During March and September, the earth will be at its maximum distance, 27 
solar radii, South and North of the sun's equator. According to Figure 4 the earth during 
these months will be at least 12 solar radii outside of the direct equatorial flow, and in 
the presumably turbulent fringe region of the coronal disk. Therefore, geomagnetic dis- 
turbances may be expected to be especially prevalent during these seasons as has been, 
in fact, observed (Reference 5). 


EFFECT OF THE SOLAR WIND ON THE LUNAR ATMOSPHERE 

The composition and extent of the lunar atmosphere are determined by the competi- 
tion between the processes of accumulation of gases from the surface and the escape 
from the atmosphere of those molecules with velocities sufficiently high to remove them 
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from the gravitational attraction of the moon. Probab e sources for the accumulation of 
surface gases are: the radioactive decay of potassium, producing argon; and residual 
volcanic activity, producing such gases as S0 2 , C0 2 , and H 2 0. Estimates of the emis- 
sion of these gases based on terrestrial data suggest a source of strength of 5 * 10 5 
cm' 2 sec' 1 for argon, and 10 10 cm' 2 sec' 1 for the volcanic gases (Reference 6). These 
estimates are, of course, highly questionable, and are used here for illustrative purposes. 

The moon has one -eightieth of the mass and one-quarter of the radius of the earth. 

Its gravitational potential is therefore one-twentieth that at the surface of the earth. As 
a result of this low gravitational potential, gases escape much more rapidly from the 
lunar atmosphere than from the earth’s atmosphere. The rate of escape also depends on 
the maximum temperature, which determines the number of atoms or molecules in the 
high-energy tail of the thermal distribution. Thermocouple measurements indicate a 
value of 370°K for this temperature (Reference 7). 


If an isothermal atmosphere and a Maxwellian distribution of particle velocities 
with an infinite mean free path are assumed, a calculation based on the balance between 
these two processes leads to the following result for the equilibrium number density of 
atmospheric molecules at the lunar surface (Reference 8): 


n 0 = jV 


rrm 

2 kT 


(• 


m M G\“ 1 m MG 
k TRj EXP kTR ' 


(24) 


In Equation 24, J is the the number of molecules ejected from the lunar 
pcy nnnnnri; T the lunar atmospheric temperature; m 
ihp mass nf the moo n: and f its radius. 


surface per 
the mass of an 


When these values of J and T are used in Equation 24, it is found that for argon, 
n 0 % 5 x 10 1 5 /cm 3 or 10' 4 atmosphere. A similar calculation for the volcanic gases 
leads to pressures greater than 1 atmosphere for SO 2 and C0 2 , and to 4 * 10 8 atmos- 
phere for H 2 0. However, the latter values will be lowered substantially when allowance 
is made for the effects of photodissociation and chemical reactions with the crust, and 
for the fact that the age of the moon (~ 4.5 x 10 9 years) does not allow enough time for 
S0 2 and C0 2 to build up to equilibrium (Reference 6). 


These estimates of the extent of the lunar atmosphere are seriously reduced by the 
effects of the solar wind. It appears, in fact, on the basis of conventional calculations of 
the probability of escape, that the solar wind will blow away all but a very small fraction 
of the lunar atmosphere. 


Consider a solar wind consisting of protons with the conventional density of 10 3 /cm 3 
and velocity of 10 8 cm/sec, or energy of 10 kev (Reference 9). In an elastic collision 
with an atom these protons will transfer an average of 1 kev of kinetic energy, which is 
sufficient for the escape of the atom. If it is assumed that every atom struck by a proton 
escapes, the rate of ejection of particles is readily calculated, and this gives for R p the 
value of the rate of ejection of particles from the atmosphere (per cubic centimeter): 

Rp = "a "p V p CT el ' (25) 
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2N 


In Equation 25, n p is the number density of protons, v p is their velocity, and a el 
is the proton-particle elastic cross section. The equilibrium particle density at the 
lunar surface is estimated by equating the rate of ejection of particles from a vertical 
column 1 cm 2 in cross-sectional area, as calculated by Equation 25, to the rate of 
injection of particles J into the column at the surface. This gives 


n 


o 


n V <y . h 

P p el 11 


( 26 ) 


where h is the scale height of the isothermal atmosphere kTR 2 /mMG. Taking a el to be 
10" 16 cm 2 , it is found that n 0 % 10 4 /cm 3 or 10" 15 atmosphere for argon. For the vol- 
canic gases, n 0 ~ 10 8 /cm 3 or 10" 11 atmosphere. This value represents an upper limit 
for these gases because photodissociation has been neglected. 


When these values are compared with the results obtained earlier, it is seen that the 
solar wind reduced the density of the lunar atmosphere by a factor of 10 11 for argon, 10 4 
for H 2 0, 10 15 for C0 2 , and 10 22 for S0 2 . 

These results depend on the assumption that the moon's magnetic field is too feeble 
to shield the atmosphere effectively from the solar wind. If an upper limit of 100 gamma 
is assumed for the seleno-magnetic field (References 10 and 11), then the magnetic pres- 
sure of such a field of this intensity will indeed have no shielding effect. 

The values quoted for argon and the volcanic gases lie to either side of an experi- 
mental estimate of 10' 13 atmosphere, which has been derived from measurements of the 
lunar ionosphere by Elsmore and Whitfield (Reference 12). 


THE MOON'S HYDROGEN ATMOSPHERE 

Gold (Reference 13) suggested that the high-velocity protons from the solar wind will 
be absorbed by the lunar surface and subsequently emitted by the surface as a gas of 
neutral hydrogen at the temperature of the point of emission. Since the thermal velocity 
of the emitted gas is considerably smaller than the incident solar wind velocity, the 
density will be correspondingly greater, in accordance with the equation of continuity, 
and in fact great enough to produce a substantial lunar atmosphere of cool monatomic 
hydrogen. 

The thermal velocity v s of the gas escaping from the subsolar point is about 
2.5 x 10 5 km/sec. This is slightly greater than the escape velocity V e = 2.4 x 10 5 km/ 
sec. It follows that a large fraction of the atoms in this atmosphere will be in the 
process of escaping. A crude estimate of the density N H may be obtained by considering 
the atmosphere to be in a state of radial flow outwards at the constant velocity V s . The 
equation of continuity gives 



and the constant of proportionality is obtained by requiring that the flux j of solar pro- 
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tons into the surface shall equal the outward flux of hydrogen atoms. Thus 



where R is the moon's radius. With J = 10 10 cm' 2 see' 1 , the quantity N H = 0.4 x 10 5 / 
cm 3 is obtained at the moon's surface. 


Summarized next is an attempt at a closer description of the hydrogen atmosphere 
taking into account the variation in the moon's surface temperature and the variation in 
the component of the solar proton flux normal to the moon's surface. 

Denote by da(P) an area element at the point P on the moon's surface. Each point 
receives particles from the incident solar stream at a rate Jj(P), and from the rest of 
the moon's surface at a rate J A (P). 

The particles absorbed from the solar stream are mostly protons. They will be 
emitted as neutral hydrogen atoms at a rate R (P) . A certain fraction F e (P) of the flux 
R(P) consists of particles that escape to infinity; and the remainder consists of parti- 
cles that return to the moon, where they are absorbed and subsequently re-emitted. 

The fact that the thermal velocities involved are comparable with the escape velocity 
implies that the average path length between emission and absorption for the returning 
particles is a sizeable fraction of the moon’s circumference. Thus, particles from a 
rather large area contribute to the rate J A (P). It therefore seems reasonable to assume 
that J A (P) is constant, J A (P) = J A . 

If it also is assumed that the atmosphere is in a steady state, two equilibrium con- 
ditions must be satisfied by the various rates. First, to a good approximation each point 
must emit as many particles as it receives: 

R(P) = J A + Ji( p )- ( 27 ) 

Second, as many particles must enter the atmosphere as leave, 

J*F e (P) R(P) da(P) = Jj^P) <MP). (28) 

where the integrals are taken over the moon's surface. 

Denote the angle between the point P and the subsolar point by 6 . Then 


Jj = J cos 0 for 0 < # < > 

Jj = 0 for 0 >-y . 


(29) 


The quantities R(P) and F e (P) can also be calculated if it is assumed that the emitted 
particles have a Boltzman distribution with density N(P). Then 
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and 


R(P) 


2 

— N(P) V(P), 

V 77 



where V(P) is the local thermal velocity, 


v<p>=y5p 


(30) 

(31) 


(32) 


and V e is the escape velocity, 

V e = VW • (33) 

If the surface temperature T(P) is known, Equations 29 to 33 convert the equilibrium 
conditions (Equations 27 and 28) into two equations for the unknowns J A and N(P). 

When the surface temperature determined by Pettit and Nicholson (Reference 14) is 
used, the curve shown in Figure 5 is obtained for 

N x 10 5 

J 

as a function of cos 8 , where 6 is the angle between P and the subsolar point S. For a 
solar wind flux J = 10 10 cm' 2 sec' 1 , there is obtained for the subsolar point 

N s = 0.39 x 10 5 /cm 3 , 

in good agreement with the more qualitative argument. On the dark side of the lunar 
surface, N D = 0.085 x 10 5 /cm 3 , with a ratio N S /N D = 4, is obtained, and also J A = 0.14 J. 



Figure 5 - The hydrogen density on 
the Lunar surface 
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THE LUNAR IONOSPHERE 

In the lunar ionosphere the dominant ionization process is charge exchange between 
solar protons and hydrogen atoms. The cross section for this process is about 10 t5 cm 
as compared with 10 ' 18 cm 2 for photoionization. In the following the ionization density 
due to charge exchange is computed upon the assumption that every ionized particle escapes 
from the moon’s atmosphere. Although some newly ionized protons may be knocked into 
gravitationally bound orbits that do not intersect the lunar surface, it is not possible to 
estimate the density of such protons without a more detailed analysis of the velocity dis- 
tribution in the hydrogen atmosphere. The ionization density computed from charge ex- 
change is therefore to be regarded as a lower limit. 

The effect on this atmosphere of charge exchange between the hydrogen atoms and 
the solar wind protons is now calculated as follows: The cross-section a for this 
process, for protons of energies between 1 and 10 kev., is between 2 x 10 and 
1 x 10- 15 cm 2 (Reference 15). The number of slow protons produced per cubic centime- 
ter per second is N H J a. Let N p denote the density of the slow protons. If the moon’s 
gravitational field and the energy transfer in the charge exchange are neglected, then it 
can be assumed that both the slow protons and hydrogen atoms are streaming radially 
outwards with constant velocity V. 

The continuity equations for these two constituents of the atmosphere are: 

72 £ <"» v > ■ - n h J" 

72 £ <»p ' Jv > * * N « J»- 

It is found that 

n h = v(t ) 2 exp --ir( r - R ) 

\ = ‘ exp -ir (r ' R) } 

Now 

J = 10 10 cm -2 sec ' 1 
V - 10 s cm/sec ; 


thus, 






X 10 - 15 

To 3 


x 


1.7 x 


10 8 


~ 1. 7 x icr 2 , 


so that for r < 100 R, it may be taken that 



1(f) 2 ■ 
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The bracketed term possesses a maximum of 1/4 at r 2R. With these values, then 




At r = 2 R , N p = 400/cm 3 and at r ^ 10 R, N p = 170/cm 3 and should just begin to be 
detectable above the interplanetary density of 100/cm 3 . This is in agreement with the 
results reported by Krassovskii (Reference 11) for Lunik III, in which the interplanetary 
ion density was found to begin to increase above the ambient value at a distance of about 

r = 10 R . 
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